A physical interpretation for the non-Hermitian Hamiltonian 
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We explore a way of finding the link between a non-Hermitian Hamiltonian and a Hermitian one. 
Based on the analysis of Bethe Ansatz solutions for a class of non-Hermitian Hamiltonians and the 
scattering problems for the corresponding Hermitian Hamiltonians. It is shown that a scattering 
state of an arbitrary Hermitian lattice embedded in a chain as the scattering center shares the same 
wave function with the corresponding non-Hermitian tight binding lattice, which consists of the 
Hermitian lattice with two additional on-site complex potentials, no matter the non-Hermitian is 
broken VT symmetry or even non-VT- An exactly solvable model is presented to demonstrate the 
main points of this article. 

PACS numbers: 03.65.Ge, 05.30.Jp, 03.65.Nk, 03.67.Bg 



I. INTRODUCTION 

In general, a non-Hermitian Hamiltonian is said to be 
physical when it can have an entirely real energy spec- 
trum. Much effort has been devoted to establish a parity- 
time (VT) symmetric quantum theory as a complex ex- 
tension of the conventional quantum mechanics [H-IH 
since the seminal discovery by Bender [![. It is found 
that non-Hermitian Hamiltonian with simultaneous VT 
symmetry has an entirely real quantum mechanical en- 
ergy spectrum and has profound theoretical and method- 
ological implications. Reseaches and findings relevent to 
the spectra of the VT symmetric systems are presented, 
such as exceptional points [§], spectral singularities for 
complex scattering potentials [10], complex crystal and 
other specific models [ll| have been investigated. At the 
same time the VT symmetry is also of great relevance to 
the technological applications based on the fact that the 
imaginary potential could be realized by complex index 
in optics |13 - [l6| . In fact, such VT optical potentials can 
be realized through a judicious inclusion of index guiding 
and gain/loss regions and the most interesting aspects as- 
sociated with VT symmetr ic sy stem are observed during 
dynamic evolution process |17H20f . 

Thus one of the ways of extracting the physical mean- 
ing of a pseudo- Hermitian Hamiltonian with a real spec- 
trum is to seek for its Hermitian counterparts [2ll - l23l ]. 
The metric-operator theory outlined in Ref. provides 
a mapping of such a pseudo-Hermitian Hamiltonian to 
an equivalent Hermitian Hamiltonian. Thus, most of the 
studies focused on the quasi-Hermitian system, or un- 
broken VT symmetric region. However, the obtained 
equivalent Hermitian Hamiltonian is usually quite com- 
plicated [1, [H| , involving long-range or nonlocal interac- 
tions, which is hardly realized in practice. 

To anticipate these problems, alternative proposals for 
the connection between a pseudo-Hermitian Hamiltonian 
and a real physics system have been suggested in the con- 
text of scattering problems [25:j . Central to that analy- 
sis was the recognition that the VT Hamiltonian may 
be used to depict the resonant scattering for an infinite 
system. It is shown that any real-energy eigenstate of 



certain VT tight-binding lattice shares the same wave 
function with a resonant transmission state of the corre- 
sponding Hermitian lattice. In such a framework, further 
questions to ask are whether the requirements of the en- 
tireness of the real eigenvalunes and the VT symmetry 
of the non-Hermitian system are really necessary. 

In this paper, we propose a physical interpretation for 
a general non-Hermitian Hamiltonian based on the con- 
figurations involving an arbitrary network coupled with 
the input and output waveguides. Relevant to our pre- 
vious discussion is the interpretation of the imagiary po- 
tentials. Based on this, we make a tentative connection 
between a non-Hermitian system and the corresponding 
large Hermitian system. It is shown that for any scatter- 
ing state of such a Hermitian system, the wavefunction 
within the center lattice always corresponds to the equal 
energy eigenfunction of the non-Hermitian Hamiltonian, 
no matter it is VT symmetric or not. Our formalism is 
generic and is not limited to the pseudo-Hermitian sys- 
tem. 

This paper is organized as follows. Section II is the 
heart of this paper which presents a formulism to reduce 
a scattering process of a Hermitionian system to the eigen 
problem of the non-Hermitian system. Section III con- 
sists of two exactly solvable examples to illustrate our 
main idea. Section IV is the summary and discussion. 



II. NON-HERMITIAN REDUCTION OF A 
HERMITIAN SYSTEM 

A typical scattering tight-binding network is con- 
structed by a scattering-center network and two semi- 
infinite chains as the input and output leads. The well- 
established Green function technique [IB-HH can be em- 
ployed to obtain the reflection and transmission coeffi- 
cients for a given incoming plane wave. The correspond- 
ing wave function within the scattering center should 
be obtained via Bethe ansatz method. In the following 
we will show that this can be done by solving a finite 
non-Hermitian Hamiltonian. In our previous work [25j . 
a VT symmetric non-Hermitian Hamiltionian has been 
connected to a physical system in the manner that any 
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real-energy eigenstate of a VT tight-binding lattice with 
on-site imaginary potentials shares the same wave func- 
tion with a resonant transmission state of the correspond- 
ing Hermitian lattice embedded in a chain. The main aim 
of this article is to answer the question of whether such 
a statement still holds for broken VT non-Hermitian or 
mm-VT lattice. In the following, we will show that a 
scattering state of the Hermitian system always has con- 
nection to the eigenstate of its non-Hermitian reduction. 
For a certain incident plane wave, the scattering problem 
of the whole infinite Hermitian system can be reduced to 
the eigen problem of a finite non-Hermitian system. 

The Hamiltonian of a typical scattering tight-binding 
network has the form 



H = H A + H B + H C 



(1) 



where 



H A = -jJ2 b t-i b i - 9AbUa A + H.c. (2) 

i=-l 
+ oo 

H B = -Jj2b\bi + i - g B b\a B + H.c. (3) 



represent the left and right waveguides and 

N 

H c = — Kija\aj + H.c. 



(4) 



describes an arbitrary TV-site network as a scattering cen- 
ter. Sites A and B are arbitrary within the network. Here 
hi, di, are boson (or fermion) operators, — Ka (we denote 
Va — —k-aa Vb = —k-bb only for the sake of simplicity) 
represents the potential at site i. Fig. [TJa) represents a 
schematic scattering configuration for an arbitrary net- 
work. 

For an incident plane wave incoming from waveguide A 
with energy E — — 2Jcos(fc), the scattering wave func- 
tion can be obtained by the Bethe ansatz method. The 
wave function has the form 



\ipk) = ^2 fibj |vac) + ^2 h i a \ l vac ) 



where 



fl = 



e ifc(I+l) +re -«(l+l) ) l g (-oo,-l] 



te ik(l~l) 



I G [l,oo) 



(5) 



(6) 



hi = hi, I G [1,N] 



Here r, t are the reflection and transmission coefficients. 
The explicit form of the Schrodinger equations for the 
waveguides Ha and Hb are 



H\i> h ) = E\^ k ) 



(7) 



admits 




(a) 




(b) 



FIG. 1. Schematic illustration of the configuration of the 
concerned network. It consists of an arbitrary graph of a 
Hermitian tight-binding network (shadow) connecting to two 
semi-infinite chains L and R as the waveguides. The wave 
function within the scattering center for a scattering state of 
the whole system is identical to an equal-energy eigen function 
of the non-Hermitian Hamiltonian which is constructed by the 
center Hermitian network with imaginary potentials added at 
the joint sites A and B. 



—Jfj-i - Jfj+i — Efj, 
0'e(-oo,-2]U[2,+oo)) 
-Jf-2 - 9 A^ A = Ef-x 
-Jfs - 9Bh B = Efi. 
From Eq. (j5]), we obtain E = —2 J cos (k) and 



(8) 



J 



h A = — (e ik +re- lk ^ 
9 A 

h B = —te- tk . 
9B 



(9) 
(10) 



Vanishing hA {hs) is beyond of our interest. From Eqs. 
© and (fTU]) , one can express the wavefunctions of two 
joints (A, B) as, 



f 1 = 9A 1 + r h A 



J e + re 



f _ 9b ikh 
h — — e n B- 



(11) 



(12) 



The explicit form of the Schrodinger equations for H c can 
be written as, 



N 



- ^ K ij h i = Eh 3' (i ^ A B) 

i 

N 

KjAhj - gAf-i = (E - Va) h A , (13) 



i^A 
N 



KiBhi - 9b fi = (E - V B ) hi 

i=£B 
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Substituting the expression for /_i and fi from Eqs. (fTTj) 
and (H2), to the above Eqs. (fT3|) . we get the following 
Schrodinger equations for the center network, 



jv 



Kijhi = Eh 3 , (j ^ A, B) 

i 
N 

^ KiAht = (E - U A ) h A , 



(14) 



N 



2J KjB^j = {E - U B ) hi 



with 



C/a = V A - 
U B = V B ~ 



l + r 



(15) 



9% ik 
J ■ 



This is equivalent to the effective non-Hermitian Hamil- 
tonian 



H = H C + (U A - V A ) n A + (U B - V B ) n B . 



(16) 



Without losing generality, we take r — \r \ e lS with \r\ < 1. 
This leads to 



lm(U A )lm(U B ) = — 



(g A g B sinkf (1 - |r| 2 ) 



J 2 



l + 2|r|cos (8 -2k) 



< 0, 



S 17) 

which means that the imaginary part of the additional 
potentials have opposite signs, one providing gain and the 
other loss. This is in accordance to the conservation law 
of the current. It is important to stress that magnitude 
of the two imaginary potentials may not equal, which 
deviates from the general understanding of an imaginary 
potential. 

The existence of the scattering solution of the Hermi- 
tian system H ensures that there must exist at least one 
real solution of T-L with eigenvalue equals to the incident 
energy E. It possesses the identical wavefunction as that 
of the scattering state within the region of the scattering 
center. Then a scattering problem is reduced to the eigen 
problem of a non-Hermitian Hamiltonian. This conclu- 
sion is an extension of our previous result [2o| . In this 
work, our formalism is generic: The central network is 
not limited to the linear geometry and the scattering is 
not restricted to be resonant transmission. Thus the scat- 
tering interpretation for the non-Hermitian Hamiltonian 
is not limited to pscudo-Hermitian system. This rigorous 
conclusion has important implications in both theoretical 
and methodological aspects. 

Likewise, if we consider the inverse scattering pro- 
cess, i.e., taking the time-reversal operation on the above 



(a) 



2N 




(b) 



FIG. 2. (Color online) Schematic illustration of the concrete 
configuration for a scattering system. A ring as the scattering 
center, connects to two semi-infinite chains L and R as waveg- 
uides with coupling —g. The on-site potentials at the connec- 
tions are Va and Vb ■ The wave function within the scattering 
center for a scattering state of the whole system is identical to 
an equal-energy eigen function of the non-Hermitian Hamil- 
tonian which is constructed by the center Hermitian ring with 
imaginary potentials U A and Ub added at the joint sites A 
and B. 



mentioned scattering process. The corresponding Bethe 
ansatz wave function has the form 



ik{i+\) 



5,ife(iH 



* p -ik(i-i) 



t*e 



n ■ 



■ l \ le (-00,-1] 

l G [l,oo) 
I G [1,M 



(18) 



with energy E = — 2 J cos k. The above conclusion still 
holds. Straightforward algebra shows that the corre- 
sponding non-Hermitian reduction is TV . In the frame- 
work of non-Hermitian quantum mechanics, IV takes an 
important role to construct a complete biorthogonal basis 
set, which has no physical correspondence. In the con- 
text of our approach, 1-0 has the same physics as H, in 
describing the scattering problem of the same Hermitian 
system. 



III. ILLUSTRATIVE EXAMPLES 

In this section, we investigate simple exactly solvable 
systems to illustrate the main idea of this article. We 
will discuss two examples which correspond to a VT 
and a non-VT non-Hermitian Hamiltonian, respectively. 
The advangtage of these examples are that the non- 
Hermiltian Hamiltonians are exactly solvable. 



A. Exactly solvable VT Hamiltonian 

To exemplify the previously mentioned analysis of re- 
lating the stationary states of a non-Hermitian VT- 
symmetric Hamiltonian to a scattering problem for a 
Hermitian one, we take the center network to be a sim- 
ple network: a uniform ring system. We start with the 
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scattering problem for a class of symmetric systems, the 
Hamiltonian can be written as 

2N 2 

H ss = -J^ a\a i+1 + H.c. - - — (m + n N+1 ) 

i—l " 

— oo +00 

-J J2 6 i-i b < - +H.c. (19) 

i=-l i=l 

— .g&^ai — .gbjaAr+i + H.c, 

where we denote the connection sites as = a± and 
as = Oiv+i- 

The corresponding non-Hermitian Hamiltonian de- 
pends on the energy E of the incident plane wave as 
well as the parameters e and g. To be concise, as an 
illustrative example, we would like to present the ex- 
actly solvable model, which are helpful to demonstrate 
our main idea. Therefore, we will focus on the following 
configurations: 

i) g 7^ V2J, E — e = e n = — 2Jcos (nn/N) where 
n £ [1,7V — 1]. Here we restrict the energy of the inci- 
dent plane wave since it will leads to the pure imaginary 
potential, thus ensures the existence of the exact solu- 
tion. Straight forward algebra shows that the problem of 
solving the Schrodinger equation is reduced to the eigen 
problem of the following non-Hermitian Hamiltonian 

27V 

n [n] = -J^ a\a i+ i + H.c. + i7„a|ai - i7 n a] v+1 ajv+i 
i=i 

(20) 

with the imaginary potential 

q 2 . TL7T . 

7« = y sin (^)' ( 21 ) 

Obviously, this Hamiltonian depicts a 2A-site ring with 
two imaginary potentials at two symmetrical sites, which 
is a 'PT-invariant Hamiltonian. Note that the magnitude 
of the imaginary potential is discrete in order to obtained 
the exact solutions. In Appendix A, it is shown that 
such lattices can be synthesized from the potential-free 
lattice by the intertwining operator technique generally 
employed in supersymmetric quantum mechanics. The 
eigen spectrum of consists of 

e j = -2 J cos (jir/N), (22) 
(je [1, N - 1] , 2-fold degeneracy) 

and two additional levels 

e± = ±v/4J 2 - ll (23) 

The eigenstates with eigenvalue Ej can be decomposed 
into two sets: bonding and antibonding, with respect 
to the spatial reflection symmetry about the axis along 
the waveguides. For the scattering problem, only the 
bonding states are involved. It shows that there always 



exists a solution in {sj} to match the energy e n of the 
incident wave. 

From Eqs. (|2"2l |2"3"|) . on can see that a pair of imagi- 
nary eigenvalues appear, i.e., the VT symmetry is broken 
when g > y/2J . In general, a non-Hermitian Hamilto- 
nian with a broken VT symmetry is unacceptable be- 
cause its complex energy eigenvalues make a hash of the 
physical interpretation. On the other hand, the VT sym- 
metry breaking was observed in optics realm experimen- 
tally [29j]. In theoretical apects, VT symmetry in non- 
Hermitian spin chain system was discussed [301 ] - From 
the point of view of this article, we note that even 
possesses a broken VT symmetry, the spectrum {ej } still 
contains the state with the energy ej — e n . It is worth 
mentioning that the broken symmetry does not contra- 
dict the interpretation of the non-Hermitian Hamiltonian 
(|2T)]) . This idencates that even the VT symmetry is bro- 
ken the non-Hermitian Hamiltonian still has physical sig- 
nificance. 

ii) g = y/2J, E = e £ [-2J, 2J]. Here we do not 
restrict the energy of the incident plane wave but the 
magnitue of g. Straight forward algebra shows that the 
problem of solving the Schrodinger equation is reduced to 
the eigen problem of the following non-Hermitian Hamil- 
tonian 

2N 

H [e] = -./y^ a\a l+x + H.c. + i^ e a\ai - i7 e aj v+1 a A r + i 
i=l 

(24) 

with the imaginary potential 

7e = ^4J 2 - e 2 . (25) 

From Appendix A, the solution of the Hamiltonian T-L^ 
has the same form of Eqs. (|22l [23]) with j n replaced by 
7 e . Here we would like to see the relation between the 
Hamiltonians and W^: Both of them come from 
the same model with different coupling constants (with 
g 7^ \[2J and g = V2J) and different incident plane 
waves (with discrete and continuous spectra). However 
they have the same structure but different values of the 
imaginary potentials. In Appendix A, we provide the 
universal solution contains that of % and H M . 

Obviously, Hamiltonian is always exact VT sym- 
metric. All the eigenvalues are real. Among them we 
can find that e± = ± -J A J 2 — 7! = ±e, one of e± equals 
to the energy of incident plane wave e and thus verifies 
the above mentioned conclusion. Furthemore, the solu- 
tion of it has the following peculiar feature: in the case 
of e = e n , i.e., the incident wave has wave vector nir/N 
(n £ [1,N — 1]), the exceptional points appear in W^. 
It is shown in Appendix A that the corresponding eigen- 
functions of e+ and e„ (en-u) coalesce. 

According to non-Hermitian quantum mechanics, in 
general, has the Hermitian counterpart which 
possesses the same spectrum. When the potential 7 e ap- 
proches 7 £n , the similarity transform that connects 
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and ff[ £ l becomes singular. The Hamiltonian be- 
comes a Jordan-block operator, which is nondiagonaliz- 
able and has fewer energy eigenstates (N — 1) than eigen- 
values (iV+1), (i.e., the lack of completeness of the energy 
eigenstates.) Such a Hamiltonian has no Hermitian coun- 
terpart [3l[. According to our analysis, one can see that 
even at the exceptional points @ the coalescing eigen- 
states still has physical significance. 

B. Exactly solvable non- VT Hamiltonian 

Now we turn to exemplify the previously mentioned 
analysis of relating the stationary states of a non- 
Hcrmitian non-7'T-symmetric Hamiltonian to a scatter- 
ing problem for a Hermitian one. We still take the center 
network as a simple network: a uniform ring system with 
uniform coupling but none on-site real potentials. The 
corresponding Hamiltonian can be written as 

2N 

#as = - J^otoi+i+H-c. (26) 
1=1 

— oo +oo 
i=-l f=l 

— Jb_-yCLi — Jb\a,N+i + H.c. 

We consider the incident plane wave with energy E — 
— 2Jcos0, where d £ (it, —it) without any restriction. 
Straight forward algebra shows that the problem of solv- 
ing the Schrodinger equation is reduced to the eigen prob- 
lem of the following non-Hermitian Hamiltonian 

27V 

H [§] = - a\a i+ i + H.c. + U A a\ai + U B a> N+1 a N+1) 

i=l 

(27) 

where the complex potentials are 



Hamiltoian. In this sense, one can conclude that a non- 
VT non-Hermitian Hamiltonian still has physical signif- 
icance. 



IV. CONCLUSION 

In summary, we have studied the connection between a 
non-Hermitian system and the corresponding large Her- 
mitian system. We propose a physical interpretation 
for a general non-Hermitian Hamiltonian based on the 
configurations involving an arbitrary network coupled 
with the input and output waveguides. We employed 
the Bethe ansatz approach to the scattering problem to 
show that for any scattering state of a Hermitian system, 
the wavefunction within the scattering center lattice al- 
ways corresponds to the equal energy eigenfunction of 
the non-Hermitian Hamiltonian. It is important to stress 
that such a physical interpretation for the non-Hermitian 
Hamiltonian is not limited to the pseudo-Hermitian sys- 
tem. As an application, we examine concrete networks 
consisting of a ring lattice as the scattering center. Exact 
solutions for such types of configurations are obtained to 
demonstrate the results. Such results are expected to be 
necessary and insightful for the physical significance of 
the non-Hermitian Hamiltonian. 
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Appendix A: Construction of 'PT-Hamiltonian by 
Interwining operator technique 

In this Appendix, we will derive the central formula 
for studying the eigen problem of the VT ring system. 



U A = 
U B = 



-J 



e MN cos [(AT - 1)0] +ismi9 
sin [(AT - 1)0] -sintf 



-Je 



2sin0, (28) 



We can see that, in general, Hamiltonian %^ is not VT 
symmetric, except in some special cases. It is hardly 
to get the analytical solution of such a Hamiltonian in 
general cases. Fortunately, what we need to do is to 
prove that the incident energy E — — 2 J cos is always 
one of the eigenvalues of the Hamiltonian. In fact, in 
single-particle basis the matrix representation J\4^ of 
the Hamiltonian (|2"T|) satisfies 



det M m +2Jcos0 = 



(29) 



according to the derivation given in Appendix B. This 
result do not depend on the pseudo-Hermiticity of the 



1. Linear Transformation 

First of all, the Hamiltonian can be decomposed into 
two independent sub-Hamiltonians 



1~l — 



N-l 



n B = -j 



N-l 
\^ fit 



i+l 



H.c. 



(Al) 



+U A a\a 1 + U B a* N+1 a N+ i, (A2) 



(A3) 



6 



Ua 



Ub 



B 



N+1 



FIG. 3. (Color online) Schematic illustration of the reduc- 
tion for a ring system by linear transformation. The top 
panel (black) represents the Hamiltonian T-L a , while the bot- 
tom panel (red) shows the Hamiltonian Hp 



with [H a , Hp] = 0, by using the following linear tranfor- 
mation: 

o?i = ai,Q!jv+i = ctjv+i, 

«j = -j= (oj + «2W+2- J ) , .? € [2, AT] , (A4) 

/3j = (oj - fl2iv+2-i) , j e [2, AT] . 

We will focus on the solution of the Hamiltonian H a . 
Typically, the solution can be obtained via Bethe ansatz 
method as shown in Ref. [28]. In this Appendix, we 
will use the intertwining operator technique to get the 
solutions in order to reveal their characteristic features. 



2. Interwining operator technique 

The intertwining operator technique is generally em- 
ployed in supersymmetric quantum mechanics, which 
provides the universal approach to creating new exactly 
solvable models. Recently, it is applied to discrete sys- 
tems in order to construct the model which supports the 
desirable spectrum [H, HI] . 

The critical idea of the intertwining operator technique 
is as the following: Consider an JV x ]V Hamiltonian Hi 
which has the form Hi = QilZx + fii, where Qi and 
72.1 represent N x (N + 1) and (N + 1) x N matices, re- 
spectively. One can construct an (N + 1) x (N + 1) new 
Hamiltonian H2 (H2 = IZiQi by interchanging the 

operators IZi and Q\. The spectrum of Hi is the same 
as that of H2 except for the energy level fi\. Iterating 
this method results in a series of Hamiltonians H 3 , H^, 
H$, ■ ■ ■ whose energy spectra differ from that of Hi ow- 
ing to the addition of the discrete energy levels {lii, ll 2 }, 

{/"l, M2, /^3}, {Ml, M2, ("3, 1*4,}, 

Our aim is to construct a "PT-invariant Hamiltonian 
i?3 by adding two energy levels E = — /i and E = fx 
(0 ^ fi ^ 2, the obtained conclusion will be extended 
beyond this region later) into the energy spectrum of a 
uniform chain system. We will show the processes of 



this construction explicitly. We start with the following 
(N - 1) x (N - 1) Hamiltonian 



JV-2 



Hi = -J2(\ n ) (n + l|+H.c. 



(A5) 



which depicts an (N — l)-site uniform chain. The spec- 
trum of Hi can be expressed as 

7Z7T 

e n = -2cos( — ),ne [1,AT- 1]. (A6) 

On the other hand, Hi can be written in the form 

Hi = QTZ - (i (A7) 

where 

N-l 



Q=^2(q n \n){n\+q n \n){n+l\) (A8) 

n=l 

N-l 

U= ^M n ) (n\+f n \n+l) (n\) 



n=l 



and 



(X = 2cosk, (k > 0) 



' n — Qn 



-in/2 



r n = q n = e lK/2 



(A9) 



Then the Hamiltonian H2 can be constructed in the form 
H 2 =TZQ-li (A10) 

N-l 

= -J2 (\n) (n + l\ + H.c.) - e" |1) (1| - e*« \N) (N\ , 

n=l 

which possesses an extra eigenvalue —2 cos k based on the 
spectrum s„. 

Next step, we repeat the above procedure based on a 
new Hamiltonian H' 2l which is obtained from H2 under 
parity operation P, i.e., 



H! 2 = P^ 1 H 2 P 

N-l 



(All) 



J2 (\n) (n + l\+ H.c.) - e- iK \1) (1| - e lK \N) (N\ 



where 



Pij — Si.N+l-j 



(A12) 



is the matix representation of mirror reflection. Note that 
H'2 and H2 have identical spectra. Accordingly, H' 2 can 
be written as the form 



H! 2 = Q'W + Li 



(A13) 



where 



7 



Q! 

K' 



N 

£ 

n=l 
N 

n=l 



(q' n \n) (n\+q' n \n) (n + l|) 
|n+l) (n|) 



(A14) 



and 



' JV 



= q' n = ie^/ 2 
= q[ = iV2e- lK ? 2 
= q' N = iV2e iK / 2 



(A15) 



where fi = 2 cosh u > 2. On the other hand, if k is 
replaced by a complex number k — tt/2 — iui, one can 
obtain a non-Hcrmitian VT symmetric Hamiltonian in 
the broken phase. In this case the two added eigenstates 
have pure imaginary eigenvalues fi = ±2isinhw. 



4. Coalescence of eigenstates 

Now we investigate the eigenfunctions in the case of 
k = k (k = nir/N,n € [1,2V — 1]). In this situation, all 
the cigenfunction can be written explicitly as 



Finally, the target Hamiltonian H3 can be constructed 
in the form 



K'Q' + ti 

N-l 



V2\l) (2| - V2\N) (N+l\ 



+2imn/s(|l)<l|-|JV + l> <JV + 1|) , 

and the energy spectrum (|2"2")l and (|2U)) and "H^can 
be obtained by adding the unit J. 



3. Eigenfunctions of H3 

Now we turn to derive the eigen functions of H3. The 
eigenfunctions of the Hamiltonians H±, H2 and H3 are 
denoted by (p n , ip n and ip n , respectively. The eigenfunc- 
tions of a uniform chain can be readily written as 



2 sin ( ^f] 
N \ N J 



(n e [l,N- 1]), (A17) 



According to the interwining operator technique of su- 
persymmetry theory, we have 



(p n (j) = K<p n ,n€ [1,2V- 



and 



Vv> {j) = n'P(p n ,ne [i,N] 

^ N+1 (j) = (-l) j e-™i 



(A18) 



(A19) 



^,j = l,2V + l 



i,ie [2,2V] 

Note that the eigenfunctions are not normalized. 

In the above, we restricted [i in the region [0, 2] for 
the purpose of obtaining H3 as a non-Hcrmitian unbro- 
ken VT symmetric Hamiltonian with imaginary poten- 
tials at the edges in the form of Eq. (|A16|) . However, 
the obtained result can be extended beyond the region. 
Actually, one can simply replace k by — iu> in all the ex- 
pressions. Then one can obtain a VT Hermitian Hamil- 
tonian with two added bound states with energy ±/i, 



ip n (1) = -n/2(-l) n sinfc 
ip n (j) = -2ie- ( - N+1 ~^ k sinfc 
ip n (2V+1) = -iV2 sinfc 



(A20) 



(A16) 
+ H.c. 



V>at-„(1) = iV2{-lf (-1)" sinfc 
l p N _ n (j) = - (-lf N -V 2t e ^+i-m sin k > (A21) 

ip N - n (2V + 1) = iy/2 sinfc 



V>2v(l)=i72(-l) n e- ifc / 2 
ip N (j) = 2ie- (N+1 -^ lk e- lk ^ 2 , (A22) 
%p N (N + 1) = iV2e- lk / 2 



*P N+1 (l) = (-iy e-^/V2 

<Pn+i (j) = e-^ _ . (A23) 

ij N+1 (N + l) = (-l) J e- ife VV2 

For odd 2V, we have ^„ oc ^jv and iptf—n °c ipN+i, which 
means the coalescence of eigenstates. Also the norms of 
the above four eigenstates vanish. For even 2V, we have 
the same conclusion except when n = N/2. In this case, 
we have tp n oc ipx — ipN-n oc tpN+i > which means the 
coalescence of the three eigenstates. Also the norms of 
the above three eigenstates vanish. 



Appendix B: Zero determinant 

In this appendix we will prove the Eq. (|29[) . Applying 
the linear transformation introduced in Appendix A, the 
22V-dimensional matrix M. ^ can be written in a diagonal 
block form, i.e., 



M W +2Jcostf 



V 
A 



(Bl) 



where T> is (2V + l)-dimcnsional, while A is TV- 
dimensional. Then we have 



det M m +2Jcos?9 = det |T>| det 



(B2) 
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Consider the the (N + l)-dimensional matrix 2?, which 
determinant D = dct 12? I has the form 



-y/2 



-V2 






-E 


-1 




-1 


-E -1 






-1 '■■ 








'•• -1 






-1 -E 






-1 



D = 



-a/2 U b -E 
(B3) 

Using cof actor expansion along the first and last rows, 
we obtain 



where Dj is the j x j determinant 



-E -1 
-1 -E -1 



D, 



1 



■• -1 
-1 -E -1 
-1 -E 



Such kinds of determinants follow the recursion formula 



(B5) 



J-2) 



(B6) 



which leads to 



D = [E 2 + U A U B -E(U a + U b ) - 4] D N ^ (B4) 
-2(U A + U B )D N _ 2 , 



Dj = 



1 _ e 2(j+l)ik 



-jik 



1 -e 



2 ;A' 



(j < N). (B7) 



Substituting the expressions for Ua, Ub, from Eq. (|28l) . 
we get that D = 0. Thus, Eq. is proved. 
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